Propagation of short pulses through a Bose-Einstein condensate 
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We study propagation of short laser pulses in a Bose-Einstein condensate taking into account 
f—^ ' dispersive effects under the conditions for electromagnetically induced transparency. We calculate 

, dispersion coefficients using typical experimental parameters of slow-light schemes in condensates. 

By numerically propagating the laser pulse, and referring to theoretical estimations, we determine 
the conditions for which dispersion starts to introduce distortions on the pulse shape. 
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Since the first demonstration of ultra-slow light propagation through a Bose-Einstein condensate (BEC)[lJ, a major 
application of slow light has been argued to be an optical memory [2]. Information storage capacity of such an optical 
memory based upon slowing laser pulses down to subluminal speeds in BECs would strongly depend on the temporal 
width of the pulses. As the temporal width of the pulses gets smaller, condensate could contain more pulses in a given 
time interval. Unfortunately, one cannot hope for using shortest pulses (such as femtosecond or even attosecond) 
available to maximize information capacity. A challenging issue to use such a short pulse in condensates is the fact 
that slow-light scheme in BEC is based upon electromagnetically induced transparency (EIT) d, 0] which occurs for 
not so broad optical frequency window. On the other hand, temporally short pulses have large frequency widths. 
Present experimental set ups use optical pulses of temporal widths in the order of microseconds. We note that EIT 
conditions can be reasonably well satisfied for pulses of widths ~ 10~ 7 which would enhance the information capacity 
, of the condensate ten times more. 

Another challenging issue in using short pulses in BECs is that their susceptibility to dispersion which may distort 
Q^ \ the pulse shape. We note that even though dispersive effects have been found small experimentally for microsecond 
CO pulses, there is no guarantee that role of dispersion would be small when we make the pulse ten times shorter. As the 
earlier theoretical models ignore the dispersion completely^, Q, it is necessary to develop a more general theory of 
slow light propagation taking into account the higher order dispersive properties of the BEC. In the present paper we 
develop such a theory, generalizing theory of dispersive EIT in a thermal uniform gasQ. Using our general theory, we 
have found that dispersive effects are insignificant for microsecond pulses used in slow-light experiments. However, 
^3 ' temporal width of microsecond is the critical width below which dispersion starts to distort the pulse shape. We show 
that at a width of 10 -7 seconds, dispersion leads to a small temporal broadening of the pulse. 

It may be noted here that study of dispersive effects in slow-short- light propagation may lead to alternative applica- 
tions such as determination of condensate temperature through measurement of broadening, pulse shape engineering, 
and frequency band narrowing. While we treat the light classically, novel effects have been predicted for slow-light 
in quantum regime, in particular a new nonlinearity regime for quantum light 0]. Our results may contribute to such 
studies as well. Furthermore, the theory in the present paper can easily be extended to situations in which BEC can 
have an enhanced nonlinear optical response of BEC that might be utilized for dispersion compensation. 

The paper is organized as follows. In Sect. 2, Electric susceptibility of an interacting BEC under EIT conditions 
is derived. In Sect. 3, wave equation including the high order dispersive terms for the EIT medium is derived. We 
present and discuss our numerical results in section 4. Finally, we conclude in section 5. 

II. EIT ELECTRIC SUSCEPTIBILITY FOR AN INTERACTING BEC 

EIT is a technique for making cancellation of induced absorption to weak probe field tuned in resonance to an atomic 
transition by applying a strong resonant electromagnetic field to couple coherently another atomic transition. Under 
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EIT refractively thick gaseous medium becomes transparent to the probe beam. In order to formulate the effect, we 
consider a three-level atom, where the atomic states are denoted by |1) , 2) and |3), A-type atom interacting with two 
lasers. Interaction of electromagnetic field with such an atom, under electric dipole approximation, can be described 
by a Hamiltonian of the form H = Hq + Hi , where 

H = 7kJi|l >< 1| + hu 2 \2 >< 2 + Kuj 3 \3 >< 3|, 

H x = -ef-E(t). (1) 

Here, H and Hi are unperturbed and interaction parts of the Hamiltonian, respectively. Probe laser, with angular 
frequency lo p , is tuned to transition 1 > to |3 > while coupling laser, with angular frequency lo c , is tuned to transition 
2 > to 1 3 >. We have assumed that only |3 > — > |1 > and |3 > — > |2 > transitions are dipole allowed. The interaction 
part of the hamiltonian for the (EIT) system can be calculated by evaluating 

Hi = -(|1)(1| + |2)(2| + |3>(3])er- £(t)(|l>(l| + |2)(2| + |3}(3|). (2) 

Denoting dipole matrix elements p 3 i — p,\ 3 = e(3|r|l) , p, 32 = p 23 = e (3| r |2} , P12 = M21 = e (l| r |2) and introducing 
slowly varying amplitudes for the electric fields E p (t) = (e/2) exp(— iojt) and E c (t) — (e/2) exp(— ioj c t), Rabi frequency 
for the coupling laser can be written as tt c = \i 32 e/%. Interaction part of the hamiltonian becomes 

Hi = -^{Mf exp(-iwt)|3>(l| 
+ ^exp(iwt)|l)(3| 

+ fj c exp(-iw c *)|3)(2| + fi c exp(iw c f)|2)(3|} (3) 

Including the decay terms phenomenologically, to take into account the finite life time of the atomic levels, to the 
Liouville Equation ihp = [H,p], density matrix equation become Q 

P = -^{H,p}~±{r, P } (4) 

where T is the relaxation matrix, which is defined by the equation < n \ T \ m >= 7„ S nm . At the same time the 
relation between T and p is given by {T, p} = Tp + pT. In terms of the matrix elements, explicit equation is given by 
il 

Pij = -T ^{H&Pkj - /?ik#kj) ~ ^ X!( rik ^ k j + ^kr k j)- (5) 

k k 

The equation of motion for the density matrix elements p^i and P21 are given by 

Phi = ~(iu3i + -^-)P3i 

- |^-exp(-icjt)(y0 33 - pu) + |Q c exp(-icj c t)|(?2i (6) 

P21 = -(«^21 + ~^)P21 

- ^^-exp(-i^)p 2 3 + ^£l c exp(iuj c t)p 3 i (7) 

The levels |3 > — > |1 > are coupled by a probe field of amplitude e and frequency to , whose dispersion and 
absorption we are interested in. The dispersion and absorption can be computed by psi- All atoms are initially in 
the ground level |1 >, p^ = I, p 2 °2 = P33 = P31 = 0- In order to solve equation ([6]), we can write the equations in 
the matrix form [9(,i? = —MR + A. where R, M and A are matrix elements. As a result R is given by R = AI~ 1 A. 

We found for the density matrix element p^i, 

ip 3 isex.p(-iu>t)(iA + T 2 /2) 

P31 = fv2~ (°) 

2^[(r 2 /2 + l A)(r 3 /2 + l A) + %] 

Finally, we determine the electric susceptibility x of BEC consisting of such three-level atoms in EIT scheme. If we 
use the relationship between electric field, atomic polarizability and (macroscopic) polarization Q, we find 

|2 



x = r ^ 01 < v " 1 z/ 1 (9) 

^ [(T 2 /2 + iA)(T 3 /2 + iA) + %' 
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where p = p{r) stands for the inhomogeneous concentration of atoms in an interacting BEC. At a given temperature 
density profile of an interacting BEC is composed of two components [T^- Density of the condensate part is evaluated 
by using Thomas-Fermi approximation, while the density of the thermal component is evaluated semi-classically in a 
harmonic trappj]]. The total ground state density is given by 

P(f> = - V(r)) &{ T c - T) + (10) 

where Uq = 4irh 2 a s /m, m is the atomic mass and a s is the s-wave scattering length, 0(.) is the Heaviside step 
function, g n (x) = TijX^ /j n , At is the thermal de Broglie wavelength, /3 = jr^f , and Tc is the critical temperature. 
We assume an external trapping potential in the form V{r) = (l/2)m(uj 2 r 2 + tu 2 z 2 ) with oj r the radial trap frequency 
and u> z the angular frequency in the z direction, p, is the chemical potential. 



III. WAVE EQUATION IN DISPERSIVE MEDIA 



In a dispersive media the relation between displacement current vector and electric field is given by 

D = e(uj)E 



(11) 



where e(u>) — cq[1 + x(u;)]. Polarization is related to the electric field through electric susceptibility via P — 6qx(u>)E. 
Susceptibility is the Fourier transformation of the response function (Green function) of the medium x(^)> x( w ) = 
J_ x(t) exp(—iujt)dt. Dispersive media are characterized by frequency dependence of susceptibility. Expanding the 
susceptibility of the dressed atom to second order about a central frequency u)q, and with P(uj — ujq) = eox( w — 
u>o)E(u — ujo), yieldsQ- 



which leads to the polarization^] 
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Maxwell equations in the medium are given by 

V • B = V D = p free 
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In addition, we have the standard constitutive relations D = eE + P, B = /ioH, J = crE. In such a medium, 
propagation of an electromagnetic wave is determined by a wave equation in the form 
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wave equation can be reduced to 
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Here, E loss term a, group velocity Vg, and dispersion coefficient hi are given by 0] 
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IV. NUMERICAL RESULTS AND DISCUSSION 



Let us first compare the relative strength of the last two terms of the wave equation. The third term describes 
the group velocity dispersion while the last (fourth) one is for the second order dispersion. We have analytically 
calculated the second order dispersion coefficient b 2 - Substituting the numerical values taken from experiment [H , where 
a s = 2.75 nm, uo r = 2ir x 69 Hz,u z = 2tt x 21 Hz, we numerically evaluated it to be in the order of b 2 ~ 10 -9 s 2 /m about 
the center of a 23 Na condensate with the total number of atoms N = 8.3 x 10 6 . For some other set of experimentally 
accessible parameters we have found b 2 ~ s 2 /m. We have also checked third order dispersion coefficient and have 
found it ~ 10~ 15 s 3 /m, which is insignificant. Taking the ratio of the third term to the fourth term, we find that the 
condition for dispersion to be significant is (t/Vg) x 10 9 <C 1, where r denotes the temporal pulse width . Subluminal 
light speeds achieved using 23 Na BECs are in the order of 1 — 100(m/s). This yields r <C 10~ 8 . This explains clearly 
laser pulses used in present slow-light experiments in condensates do not suffer from dispersion. They are however not 
the optimal choice for a higher capacity optical memory. One could use a ten times shorter pulse so that r ~ 10 -7 . 
Such a pulse would be influenced by little dispersion in the condensate and would be an optimal choice to use in 
optical memories. We present our numerical simulation results in Figs[l][2] for the propagation of pulses of widths 
10 _6 s and 10 _7 s, respectively. In the figures we use dimensionless scaled position (z) and time variables. Position is 
scaled by the condensate radius while time is scaled by the pulse width. Diameter of the BEC is ~ 50 /im. Effect of 
second order dispersion is a small broadening in the temporal pulse width as demonstrated in Fig|2l where, behavior 
of the pulse about the center of the condensate, where the dispersive effects are most influential, is shown. 

Detailed investigation of behavior of even shorter pulses require more intensive numerical studies and will be 
published elsewhere. Wave equation is solved via finite difference Crank-Nicholson space marching scheme. The 
Crank-Nicholson scheme is less stable but more accurate than the fully implicit method; it takes the average between 
the implicit and the explicit schemes [Hj]. We use forward difference scheme for the position in order to calculate the 
next step and we apply central difference scheme for the time. We discretized the wave equation and a set of linear 
equations are solved at each step to find E t+1 where i represents the position z. At each step we find E t+1 but this 
method is too slow. Accordingly, we use Thomas algorithm [l2j which is equivalent than gaussian elimination method 
for tri-diagonal matrices. 
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FIG. 1: Contour plot of the propagation of an EIT 
pulse of width r = 10" 6 s through a 23 Na BEC. Time 
and position are in dimensionless units as explained 
in the text. 
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FIG. 2: Same with FigQ]but for r = 10~ 7 s. 



V. CONCLUSION 



Ultra-slow light, which is achieved via EIT in BECs, has an effectively long optical path lengths within the con- 
densates. Therefore dispersive effects, in particular temporal broadening of the pulse, can be seen within the short 
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(micrometer) length of BEC. Normally, one would need much longer distances to be covered with the pulse propa- 
gating in other more usual medium such as an optical fiber. Though dispersion has no effect on the slowing down the 
pulse, slow group velocity makes the dispersion visible in such a short distance. 

We examined propagation of short laser pulses in a Bose-Einstein condensate taking into account the dispersive 
effects under the conditions for electromagnetically induced transparency. Our calculation of the high order dispersion 
coefficients using typical experimental parameters of slow-light schemes show that dispersive effects start to become 
influential for pulses whose widths are about 10~ 7 s. Present experiments uses microsecond pulses which are below 
that dispersion limit. 
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